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QPT1CAL WAVEGUIDE FILTERS 
Priority is being claimed from a U.S. Provisional application filed June 22, 2001 
entitled "Optical Waveguide Filters", agent's docket No, AP875USP. 

This is a Continuation-in-Part of U.S. patent applications No, 09/742,422 filed 
5 December 22, 2000 and 09/629,816 filed July 31, 2000. 
FIELD OF THE INVENTION 

The present invention relates to optical devices and especially to optical waveguide 
filters for integrated optics. 

|0 BACKGROUND OF THE INVENTION 

This specification refers to several published articles. For convenience, the articles are 
cited in full in a numbered list at the end of the description and cited by number in the 
specification itself. The contents of these articles are incorporated herein by reference and the 
reader is directed to them for reference. 
|5 In the context of this patent specification, the term "optical radiation" embraces 

electromagnetic waves having wavelengths in the infrared, visible and ultraviolet ranges. 

The terms "finite' 1 and "infinite" as used herein are used by persons skilled in this art to 
distinguish between waveguides having "finite 1 ' widths in which the actual width is significant 
to the performance of the waveguide and the physics governing its operation and so-called 
20 "infinite' 1 waveguides where the width is so great that it has no significant effect upon the 
performance and physics or operation. 

The theory of optical filters is well known in the art. Filters often are an integral part 
of optical systems. With the growing demand for dense-wavelength-division-multiplexing 
(DWDM) systems, the need for better filters with improved spectral characteristics (low 
25 sidelobe levels, narrow bandwidths, and high reflectance) is increasing. The ideal reflection 
filter response is characterized by unity reflection in the reflection band, fast roll-off of the 
response in the transition region and zero transmission outside of the reflection band. All 
existing filter technologies have non-ideal characteristics. Close to unity reflection is typically 
achieved within the band of interest, while out-of-band regions suffer from high sidelobe 
30 levels. High sidelobe levels and non-unity reflection in the reflection bandwidth can cause 
inter-channel cross talk in closely spaced DWDM channels. To avoid this effect, channel 
spacing is increased, reducing the capacity of communication networks. To remedy this issue, 
methods to improve filter characteristics are currently being sought. 
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Many filter technologies currently exist in the art, including multi-layer interference 
filters or mirrors, fiber Bragg gratings, metallic gratings, and corrugated-waveguide gratings. 
Generally, these filter-types employ Bragg reflection. None of these technologies are without 
limitations that impede on the ideality of the spectral characteristics of their filters. Basic filter 

5 design concepts are found throughout the an and apply to all of the aforementioned 

technologies. The concept of Bragg reflection is linked to the period or pitch p of the grating 
structure, which dictates the particular wavelength, or Bragg wavelength Xg, to be reflected. 
Reflection within the structure is created by the constructive addition of micro-reflections 

^ occurring at each interface, as caused by a refractive index perturbation introduced along the 

# length of the structure. A main difference among the existing technologies in the art is the 

y ; physical construction of this refractive index perturbation. 

f I The multi-layer filter or mirror is composed of a superposition of thin films with 

y! different refractive indices and layer thicknesses /. The proper selection of the thickness of 
^ each film arranged in a periodic, chirped or other fashion would create a strong reflection at a 
3tj particular wavelength and angle of incidence. 

fif Fiber Bragg gratings (FBG's) are created by photo-inducing a variation or perturbation 

|Jj of the refractive index in the core of an optical fiber constructed from appropriate materials. 
F * ? Numerous fabrication methods exist; the foremost method is based on the concept of a phase 

mask, A phase mask is a silica plate onto which a corrugated grating structure is etched. The 
20 design of the phase mask pitch is directly related to the Bragg period of the grating that will be 

photo-induced into the fiber core. A large body of work exists on the design and construction 

of FBG's. Despite the apparent maturity of this technology, fiber Bragg gratings are limited 

by fabrication issues and non-ideal out-of-band characteristics. 

Corrugated-waveguide gratings are another filter technology. They have found great 
25 use in the field of distributed-feedback (DFB) lasers. A corrugated grating, with pitch p and 

variable profile (square, sinusoidal, or other) is etched into a semiconductor waveguide layer, 

patterned using photolithography techniques. 

A corrugated metal grating comprises such a semiconductor grating coated by a thin 

metal film (several nm thick) to form a coated corrugated-waveguide grating or metal grating. 
30 The art of metallic gratings has emerged from the discovery of absorption anomalies, which 

are caused by the excitation of surface-plasmon modes at a single semiconductor-metal 

interface. The patterning of the metal layer in the form of a metal grating improves the 

coupling and mode selection of an external electromagnetic excitation to a surface-plasmon 
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mode. These devices axe also being investigated for their photonic band-gap structures. 

Reference [1] discloses a Quantum-Cascade (QC) laser which employs a metal grating 
and makes use of surface-plasmon modes excited at a single "infinitely wide" metal- 
semiconductor interface. The grating is included within the semiconductor laser structure to 
5 enhance mode selection, that is, create single mode selection through Bragg reflection, The 
grating comprises deposited strips of titanium covered by a thick evaporated layer of gold to 
create a metal-grating structure with alternating stripes of pure Au and Ti/Au defining the 
refractive index perturbation. 

At present, none of the above-described technologies are without limitations. 

I 

jO SUMMARY OF THE INVENTION 

if? An object of the present invention is to provide an alternative to the above-mentioned 

m filter design technologies, and/or to mitigate some of the limitations of the prior art. 

According to one aspect of the present invention, there is provided a grating, suitable 

© for filtering optical radiation, comprising a plurality of concatenated grating sections, physical 

HI characteristics of each section differing from physical characteristics of each adjacent section 
thereby defining a transition therebetween, at least some of the sections each comprising a 

U waveguide structure formed by a thin strip (100) of a material having a relatively high free 
charge carrier density surrounded by material having a relatively low free charge carrier 

20 density, the strip having finite width (W) and thickness (t) with dimensions such that optical 
radiation having a wavelength in a predetermined range couples to the strip and propagates 
along the length of the strip as a plasmon-polariton wave, said wave being partially reflected at 
the transition between said waveguide structure and the following said adjacent section, the 
arrangement being such that reflections at the different said transitions along said grating add 

25 constructively. 

The $trip may comprise for example, a metal or a highly doped semiconductor The 
surrounding material may comprise, for example, an insulator or undoped or lightly doped 
semiconductor, 

A plasmon-polariton filter or grating embodying this invention may be constructed by 
30 patterning a section of the waveguide strip, that is, varying its width w along the direction of 
propagation to create a physical perturbation in the waveguide over a certain section of length 
L. Another grating embodiment introduces a pattern of narrow metal gaps of appropriate size 
between metal strips over a certain section of the strip of length L to create the physical 
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perturbation. The 'gaps' may be filled by either the surrounding material, another material, e.g. 

dielectric, or another strip having a different high free charge carrier density. The pattern may 

take any form that adheres to the constraints of the applied fabrication method. 

According to a second aspect of the invention, there is provided a grating comprising a 
5 series of cells each cell comprising two grating sections, said series comprising a first set of 

cells (A t , A* ... A s ) and a second set of cells (A/, A 2 ' 3 ... As'), the two sets of cells being 

different from each other and interleaved alternately cell by cell. 

In a preferred embodiment of this second aspect of the invention, the first set of cells is 

equivalent to the second set of cells transposed longitudinally, 
IB According to a third aspect of the invention, there is provided a method of designing a 

yg grating suitable for filtering optical radiation within a specified range of wavelengths and 
l n formed from a waveguide strip surrounded by a dielectric material, the method comprising the 
W steps of 

tij (i) using a numerical analysis method, deriving for said specified wavelengths, a waveguide 
fc strip of a particular material, and a particular surrounding dielectric material, normalized phase 
J! constant ((3/p 0 ) and normalized attenuation constant (a/(J 0 ) for a particular waveguide strip 
M thickness and each of several waveguide widths, or for a particular waveguide width and for 
2 each of a plurality of waveguide thicknesses; 

(ii) determining a particular structure for the grating as comprising a series of strips having a 
20 predetermined overall length, adjacent strips in the series having different widths, or a series of 

strips all having the same width and with spaces between adjacent ones of the strips, or a 
series of strips having spaces between adjacent strips, alternate strips having different widths, 
and selecting for each of said strips a particular length; 

(iii) using the normalized phase constants and normalized attenuation constants derived in step 
25 (i), obtaining the complex effective refractive index (f^ = p/p 0 - ja/%) for each of said strips in 

said series; 

(iv) constructing an equivalent stack of dielectric slabs, each slab taking on the complex 
effective refractive index of the corresponding strip in said series of strips; and 

(v) deriving the optical response of the equivalent stack. 

30 The optical response may be derived using a transfer matrix method or coupled mode 

theory. Where the grating is uniform, however, the optical response also may be derived using 
the Bloch theorem. 

Thus, the spectral behavior of the plasmon-polariton gratings may be determined 
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through a foil-wave Transfer Matrix Method algorithm that incorporates the results of a 
Method of Lines (MoL) modeling tool that solves for the complex refractive index of the 
fundamental mode supported by the finite width symmetrical meta) waveguide sections. The 
plasmon-polariton gratings are modeled as an equivalent stack of thin dielectric slices. The 
5 group delay of the grating may be found numerically by applying a finite-difference formula of 
0(h 4 ), 

Numerous design architectures are supported by the plasmon-polariton gratings, 
among them: uniform periodic gratings; non-periodic or chirped gratings; step-chirped 
gratings; interleaved non-periodic or chirped gratings; interleaved periodic gratings; apodxzed 
§0 structures and higher order structures. All of the above mentioned design architectures are 
disclosed in full hereafter in the detailed description of the present invention 

Other advantages and features of the present invention will be readily apparent from 
the following detailed description, examples of preferred embodiments, the drawings and 
claims. 

BRIEF DESCRIPTION OF THE DRAWING 

Figure 1(a) is a cross-sectional illustration of a symmetric waveguide structure in which the 
core is comprised of a metal film of thickness /, width w and permittivity embedded in a 
cladding or background comprising an "infinite" homogeneous dielectric having a permittivity 

20 

Figure 1(b) shows the same structure in plan view; 

Figures 2(a) and 2(b) illustrate dispersion characteristics with thickness of the fundamental 
symmetric mode ss% for different waveguide widths w. (a) Normalized phase constant; (b) 

Normalized attenuation constant; 
25 Figure 3 is a plan view of a grating formed by a series of cells each comprising two waveguide 
sections having different widths and lengths; 

Figures 4(a) and 4(b) are pictorial views of typical gratings embodying this invention. 
Figure 5 illustrates an equivalent dielectric stack model used in simulating gratings embodying 
the present invention. 

30 Figures 6(a), and 6(b) illustrate the effects of the duty cycle on the spectral characteristics of a 
uniform periodic grating. The duty cycle is the fraction of the metal strip length to the period 
length. Figure 6(a) is a representation of the intensity reflectance with respect to N, where N 
is the number of periods, for various values of duty cycle in the range 25% to 75%. Figure 
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6(b) shows a 3-D representation of the variation of the Intensity Reflection with respect to N 
for various values of duty cycle. Each curve of Figure 6(b) is for one duty cycle value where 
the increment was set to 5%, 

Figure 7 illustrates the variation of the magnitude of the Intensity Reflection with N, the 
5 number of periods, for uniform periodic embodiments of the present invention, 

Figure 8 is a plan view of a periodic waveguide structure formed by a series of cells each 
comprising two opposed trapezoidal waveguide sections; 

Figure 9(a) is a cross-sectional illustration of a square waveguide structure in which the core is 
o comprised of a metal film of thickness width where t*w, and permittivity % embedded in 
'% a cladding or background comprising an "infinite" homogeneous dielectric having a 

permittivity z x ; 

hi Figure 9(b) is a plan view of a periodic waveguide structure formed by a series of cells each 
£ | comprising two waveguide sections having different widths and lengths where the waveguide 

sections have a square cross-section; 
3s Figures 10(a) and 10(b) illustrate a cubic-spline interpolation of the dispersion characteristics 

with waveguide width for a rectangular cross-section waveguide; Figure 10(a) Normalized 
O phase constant; Figure 10(b) Normalized attenuation constant; 

Figure 1 1 is a plan view of a uniform periodic grating formed by a series of cells each 

comprising two sections; a metal strip and a gap; 
20 Figures 12(a), 12(b) and 12(c) illustrate the spectral response and group delay for a uniform 

periodic grating of the present invention. Figure 12(a) illustrates the intensity reflectance and 

reflectance on a linear scale; Figure 12(b) illustrates the intensity reflectance on a logarithmic 

scale; Figure 12(c) illustrates the group delay on a linear scale; 

Figures 13(a) and 13(b) are plan views of an apodized-uniform and unapodized periodic 
25 grating respectively formed by a series of cells each comprising two waveguide sections 

having different widths and lengths; Figure 13(a) illustrates the sine apodized uniform periodic 

grating Figure 13(b) illustrates the unapodized uniform periodic grating; 

Figures 14(a) and 14(b) illustrate the sine apodized widths and effective refractive index 

respectively of a uniform periodic grating formed by a series of cells each comprising two 
30 waveguide sections; Figure 14(a) illustrates the width variation of the waveguide sections 

along the grating length; Figure 14(b) illustrates the desired sine profile of the effective 

refractive index along the grating length; 

Figures 15(a), 15(b) and 15(c) illustrate the spectral response and group delay for a uniform 
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sine-apodized periodic grating of the present invention; Figure 15(a) illustrates the intensity 
reflectance and reflectance on a linear scale; Figure 15(b) illustrates the intensity reflectance on 
a logarithmic scale; Figure 15(c) illustrates the group delay on a linear scale; 
Figure 16 is a plan view of a uniform interleaved periodic grating formed by a series of super 
5 cells each comprising four sections of equal length; two waveguide sections having different 
strip widths and two gaps; 

Figures 17(a), 17(b) and 17(c) illustrate the spectral response and group delay for a uniform 
interleaved periodic grating of the present invention; Figure 17(a) illustrates the intensity 
reflectance and reflectance on a linear scale; Figure 17(b) illustrates the intensity reflectance on 
|§ a logarithmic scale; Figure 17(c) illustrates the group delay on a linear scale; 

Figure 18 is a plan view of a linearly chirped grating formed by a series of cells each 
[J comprising a waveguide section and a gap. The period or length of the sections is varied 
Ljj linearly between adjacent cells and each cell supports a certain duty cycle; 
Z I Figures 19(a), 19(b) and 19(c) illustrate the spectral response and group delay for a linearly 
|5 chirped grating of the present invention; Figure 19(a) illustrates the intensity reflectance and 
Sj reflectance on a linear scale; Figure 19(b) illustrates the intensity reflectance on a logarithmic 

scale; Figure 19(c) illustrates the group delay on a linear scale; 
W Figure 20 is a plan view of a linearly step-chirped grating formed by a series of M uniform 

periodic sections of N cells each comprising a waveguide section and a gap where the lengths 
20 of these sections are varied linearly between adjacent periodic sections; 

Figures 21(a), 21(b) and 21(c) illustrate the spectral response and group delay for a linearly 
step-chirped grating of the present invention; Figure 21(a) illustrates the intensity reflectance 
and reflectance on a linear scale; Figure 21(b) illustrates the intensity reflectance on a 
logarithmic scale; Figure 21(c) illustrates the group delay on a linear scale; 
25 Figure 22 illustrates a chirped interleaved grating; 

Figures 23(a) and 23(b) are plan views of chirped gratings formed by a series of cells each 
comprising one waveguide section and a gap of different lengths which, when interleaved, 
would form the grating of Figure 22. 

Figures 24(a), 24(b) and 24(c) illustrate the spectral response and group delay for a linearly 
30 chirped interleaved grating of the present invention; Figure 24(a) illustrates the intensity 

reflectance and reflectance on a linear scale; Figure 24(b) illustrates the intensity reflectance on 

a logarithmic scale; Figure 24(c) illustrates the group delay on a linear scale; and 

Figure 25 is a block diagram depicting the design sequence for gratings embodying the present 
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invention. 

DETAILED DESCRIPTION 

Plasmon-polariton filters or gratings embodying this invention are useful for the 
5 processing and filtering of light, especially at communications wavelengths near 1550 nm. 

Devices of this invention are based on the waveguide configuration shown in Figures 1(a) and 
1(b) and disclosed in co-pending Internationa! patent application PCT/CA 00/01525 and in 
references [2] through [7] to which the reader is directed for reference and which are 
incorporated herein by reference. 
10 The design procedure for a grating using such a waveguide structure will be described 

first, followed by specific examples of such gratings. The design procedure is depicted in 
Figure 25 and begins with an analysis of the waveguide structure, using Method of Lines, to 
determine dispersion characteristics. 

It should be noted that the invention is not limited to the use of the Method of Lines, 
15 however. Other numerical analysis methods may be used instead; for example, the Finite 
Element Method, the Finite Difference Method or any other suitable method. 

METHODOLOGY 

Method of Lines (MoL) Numerical Technique 

20 The symmetric waveguide structure used to construct the embodiments of this 

invention, shown in Figure 1, comprises a lossy metal film of thickness t, width u> and 
equivalent permittivity % surrounded by a cladding or background comprising an "infinite" 
homogeneous dielectric of permittivity z x . The Cartesian coordinate axes used for the 
analysis are also shown with propagation taking place along the z axis, which is out of the 

25 page. Suitable materials for the waveguide strip include, but are not limited to, gold, silver, 
copper, aluminum, and highly n- or p-doped GaAs, inP or Si, while suitable materials for the 
surrounding material include, but are not limited to, glass, quartz, polymer and undoped or 
very lightly doped GaAs, InP or Si. Particularly suitable combinations of materials include Au 
for the strip and Si0 2 for the surrounding material. 

30 It is assumed that the metal region shown in Figure 1(a) can be modelled as an electron 

gas over the wavelengths of interest. According to classical or Drude electron theory, the 
complex relative permittivity of the metal region is given by the well-known plasma frequency 
dispersion relation [8]: 
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0) 



where co is the excitation frequency, © p is the electron plasma frequency and v is the effective 
electron collision frequency, often expressed as v = 1/t with x defined as the relaxation time of 
electrons in the metal. When a> 2 + v 2 < co] (which is the case for many metals at optical 

wavelengths) a negative value for the real part e r 2 is obtained, implying that plasmon-polariton 
modes can be supported at interfaces with normal dielectrics. 

The modes supported by the structures are obtained by solving a suitably defined 
boundary value problem based on Maxwell's equations written in the frequency domain for a 
lossy inhomogeneous isotropic medium. Uncoupling Maxwell's equations yields the following 
If) time-harmonic vectorial wave equations for the E and H fields: 
VxVxE - co 2 £(x,y)nE = 0 (2) 
Vx( e (x,y)' ) VxH)-o>>H = 0 (3) 

where the permittivity £ is a complex function of cross-sectional space, and describes the 
waveguide structure. For the waveguide structures disclosed herein, }x is homogeneous and 
i taken as the permeability of free space n 0 . 

Due to the nature of the numerical method used to solve the boundary value problem, 
the implicit y dependence of the permittivity can be immediately removed since any 
inhomogeneity along y is treated by dividing the structure into a number of layers that are 
homogeneous along this direction, and suitable boundary conditions are applied between them. 
20 The two vectorial wave equations (2) and (3) are expanded in each layer into scalar 

wave equations, some being coupled by virtue of the remaining inhomogeneity in e along x. 
Since the waveguide structure under consideration is invariant along the propagation axis 
(taken to be in the +z direction), the mode fields vary along this dimension according to € yz 
where y = a + j(3 is the complex propagation constant of the mode, a being its attenuation 
25 constant and J3 its phase constant, Substituting this field dependency into the scalar wave 
equations, and writing them for TE* (E, = 0) and TM* (H x - 0) modes while making use of 
V-[e(x)£] = 0 and V-H = 0 accordingly, yields simplified and uncoupled scalar wave equations 
that are readily solved. The Ey component of the TE* modes must satisfy the Helmholtz wave 



equation: 



30 




(4) 
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and the Hy component of the TNf modes must satisfy the Sturm-Liouville wave equation; 



s(x)ac" y 



+ ^?H™ + [r 2 + orns{xj\H™ = 0 (5) 



The superposition of the TE* and TM" mode families then describes any mode 
propagating in the structure analyzed. The electric and magnetic field components resulting 
from this superposition are given by the following equations: 
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In order to obtain a mode of propagation supported by a waveguiding structure, the 
Helmholtz and Sturm-Liouville wave equations (4) and (5), along with the field equations (6)- 
(1 1), must be solved for the propagation constant y using appropriate boundary conditions 
15 applied between layers and at the horizontal and vertical limits. 

The power confinement factor is defined as the ratio of mode complex power carried 
through a portion of a waveguide's cross-section with respect to the mode complex power 
carried through the entire waveguide cross-section. Formally it is expressed as: 



cf = 



(12) 



20 



where A, is usually taken as the area of the waveguide core and A„ implies integration over the 
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entire waveguide cross-section (which can be all cross-sectional space for an open structure) 
or the entire cross-sectional computational domain. S z refers to the z component of the 
Poynting vector: 

S. -- \{EJT, - E y Hl) (13) 

i and H * denotes the complex conjugate of H xr The spatial distribution of a component of 

the Poynting vector is easily computed from the spatial distribution of the relevant electric and 
magnetic mode field components. 

The boundary value problem governed by equations (4) to (1 1) is solved by applying 
the Method of Lines (MoL). The MoL is a well-known numerical technique and its 
lb application to various electromagnetic problems, including optical waveguiding, is well- 
established [9]. The MoL is rigorous, accurate and flexible. It can handle a wide variety of 
waveguide geometries, including the structures at hand. The method is not known to generate 
spurious or non-physical modes. The MoL formulation used herein is based on the 
formulation reported in [10], but simplified for isotropic media, as prescribed by equations (4) 
js - (1 1) and reported in [ 1 1]. Except for a 1-D spatial discretization, the method is exact. 

The main idea behind the MoL is that the differential field equations governing a 
waveguiding problem are discretized only as far as necessary so that generalized analytic 
solutions can be applied to derive a homogeneous matrix problem describing all modes 
supported by the structure. This approach renders the method accurate and computationally 
20 efficient since only N-l dimensions must be discretized to solve an N dimension problem. In 
the case of a two-dimensional (2-D) waveguiding structure, this means that only one spatial 
dimension needs to be discretized. The main features of this procedure, as applied to a modal 
analysis problem, are described below. 

The x axis and the function s(x) are discretized using two shifted non-equidistant line 
25 systems, parallel to the^ axis. 

The differential operators d/dx and d 1 ld 1 x in the wave and field equations are replaced 
by finite difference approximations that include the lateral boundary conditions. 
The discretized wave equations are diagonalized using appropriate transformations 
matrices. 

30 • The diagonalization procedure yields in the transform domain two systems of 

uncoupled one-dimensional (1-D) differential equations along the remaining dimension 
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(in this case along the y axis). 

These differential equations are solved analytically and tangential field matching 



5 



trr \ 



20 



25 



30 



conditions are applied at interfaces between layers along with the top and bottom 
boundary conditions. 

The last field matching condition, applied near the center of the structure, yields a 
homogeneous matrix equation of the form G(y)e = 0 which operates on transformed 

tangential fields. 

The complex propagation constant y of modes is then obtained by searching for values 



Once the propagation constant of a mode has been determined, the spatial distribution 

of all six field components of the mode are easily generated. 
• A mode power confinement factor can be computed by first computing the spatial 

distribution of S z which is then integrated according to Equation (12) 

The open structure shown in Figure 1(a) is discretized along the x axis and the 
generalized analytic solution applied along thty axis. The physical symmetry of the structures 
is exploited to increase the accuracy of the results and to reduce the numerical effort required 
to generate the mode solutions. For the symmetric structure shown in Figure 1(a), this is 
achieved by placing either electric wall (E^ = 0) or magnetic wall (H^ = 0) boundary 
conditions along the x and y axes. The remaining horizontal boundary conditions are placed at 
infinity and the remaining lateral boundary condition is either placed far enough from the guide 
to have a negligible effect on the mode calculation, or a lateral absorbing boundary condition 
is used to simulate infinite space, depending on the level of confinement observed in the 
resulting mode. 

Mode Characteristics 

The modes supported by the symmetric waveguide strip are termed plasmon-polariton 
modes, for which, the fundamental mode is the ss% mode, see reference [5]. The dispersion 

characteristics with strip thickness and width of the fundamental symmetric ss% mode are 

shown in Figures 2(a) and 2(b) for a rectangular Au strip embedded in infinite Si0 2 (£ f AVJ = - 
13 1.9475-312.65, s rSi02 = 2.085136) at a free-space wavelength of \ * 1550 ran. Figure 2(a) 
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shows the normalized phase constant and Figure 2(b), the normalized attenuation constant. 
The ssl mode possesses all six field components where the main transverse electric field 

component is the Ey component as long as the width of the waveguide is greater than the 
thickness. For a wide and sufficiently thin film (W-4^m, t~20nm), the main field components 
5 associated with this mode are the Ey and H x fields, thus the mode resembles a TEM 

(Transverse ElectroMagnetic) wave propagating along the film length. This mode is thus 
TEM-like in nature and not TM (Transverse Magnetic) as is the case for an interface between 
semi-infinite materials having positive and negative dielectric constants. 

0) Transfer Matrix Method (TMM) 

f!! Figure 3 is a schematic plan view of a plasmon-polariton filter or grating embodying 

UJ this invention. Figure 4(a) illustrates such a grating constructed by patterning a section of the 
f 3 1 metal waveguide strip, that is, varying its width, to create a perturbation of the mode effective 

refractive index over a certain length. The grating shown in Figure 4(a) can also be 
% S constructed using two different metal strips, where the wider sections are Au and the narrow 
!y sections are Ag, for example. 

f 1 Figure 4(b) illustrates a second embodiment in which a series of metal strips of uniform 

width are spaced apart by narrow gaps of appropriate length along a certain section of 
waveguide to create a perturbation of the mode effective refractive index. The gaps are filled 

20 by either the surrounding dielectric material or another dielectric. In essence, from a 

theoretical perspective, in the grating of Figure 4(b), the gaps constitute elements having zero 
width. The variation of the mode effective refractive index p/p 0 with respect to strip width for 
a particular strip thickness is clearly seen in the dispersion curves of Figure 2(a) As the strip 
width increases from 1 urn to 8 \im, the attenuation and phase constants increase. The 

25 combination of two or more sections of strips with different widths, as in Figure 4(a), or a 
combination of a metal strip with a short gap filled by a dielectric material, as in Figure 4(b), 
will create a perturbation of the mode effective refractive index along the structure length. 
The patterned perturbations may take various forms, resulting in varying functionality, as will 
become apparent from the description of design architectures hereinafter. 

30 The reflectance of either of the gratings of Figures 4(a) and 4(b), or other gratings to 

be described hereafter, was found by applying a full-wave Transfer Matrix Method (TMM) 
approach. The formulation is based on Maxwell's equations and a detailed description of the 
approach is presented in reference [12]. It was preferred, in order to adapt the TMM 
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approach, to model the embodiments of the present invention as an equivalent stack of thin 
dielectric slices, each slice corresponding to one of the grating elements, i.e. either a metal 
strip or a gap, as shown in Figure 5. Each slice takes on the complex mode effective refractive 
index, associated with the fundamental mode propagating in the corresponding waveguide 
5 section of the grating. The complex effective refractive index of the fundamental mode 
supported by the waveguide is obtained using the MoL as detailed above. The complex 

effective index was defined as n m = ^ - j ^ where fl/p 0 and a/p 0 are the normalized 

phase and normalized attenuation constants, respectively, as obtained using the MoL and 
3 plotted in Figures 2(a) and 2(b) as a function of the strip width and thickness. The equivalent 
TO stack model adopted herein supports many modes, either a TE (Transverse Electric), a TM 
(Transverse Magnetic) or a TEM (Transverse EiectroMagnetic) mode. However, the 

fundamental ss^ mode supported by the rectangular cross-section waveguide is TEM-like in 

nature as described above. A TEM mode is excited in the equivalent stack model when a 
plane wave is normally incident upon the stack along the +x-axis as shown in Figure 5. 
|s The H-field and E-field that satisfy Maxwell's equations in each layer m have the form, 

H= yH ym (x)e J " ( (14) 

E = xE Xm (x)e Je>t + zE 7m (x)e J * (15) 

where m = 0, 1, . . , N,s. Using Maxwell's curl equations written in phasor form, we proceed 
to find the field distribution of the field components in each layer. We have, 
20 V x £=-; Wr W w (16) 

VxH=j<DS 0 s r J m (17) 

Expanding the curl, we can solve for the field components in each layer using the 
general solutions of equations (14) and (15). Once the field components are known, we can 
relate the field amplitudes over each interface throughout the structure under consideration; 
25 this is achieved by ensuring that the tangential field components are continuous across each 
interface. The boundary conditions for the tangential field components between two adjacent 
slices are, 

>ix(£ m+l -£„)=0 (18) 
»x(# w+i -tf m )=0 (19) 
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Applying the tangential field boundary conditions at each interface in the dielectric 
stack allows us to relate the amplitude of the electric field at the input to the amplitude of the 
electric field at the output of the stack. From the application of the boundary conditions, the 
electric field amplitudes are related by 



e: 



= M 



El 



(20) 



where M is a 2x2 matrix given by 



M = 



M 2l M 22 



n aw 



O, (21) 



where N is the number of layers in the stack, D are called dynamical matrices for the TM wave 

1 \ 



in each layer and are given by D, = 



for normal incidence, P is the propagation 



matrix through the bulk of a slice of the stack and is given by P = 



0 e 



0 



where d is 



the length of the slice and - — n t is the propagation constant along the x-axis in each 
slice. 

The TMM approach disclosed above sequentially calculates the amplitude of the 
reflected E-field at each interface, between two different dielectric slices, throughout the entire 
equivalent dielectric stack model of the grating to obtain the input reflectance. The 
reflectance, r, is defined as the ratio of the magnitude of the reflected to incident E-field, and is 
given by, 

where R is the intensity reflectance and <p is the phase shift of the grating structure. 

An alternative formulation for the special case of uniform periodic grating structures is 
based on a Bloch-Wave formulation. The Bloch theorem exploits the periodicity of the 
structure to simplify the formulation through the introduction of a Bloch wave number, K 
For a detailed description of the application of the Bloch theorem see reference [12], The 
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theorem allows the introduction of an eigenvalue solution to simplify the formulation of 
relating the amplitude of the electric field at the input of the stack to the amplitude of the 
electric field at the output of the stack model In this situation, the reflectance is given by, 




where the unimodular matrix that relates the input and output electric field amplitudes is 
defined as, 





'A B' 










C D 





AU N _, - U N . 



BU 



cu 



N-l 



(24) 



111 



where U N = \ / 

smKA 



The group delay is another important performance characteristic of gratings. The 
group delay t of an input signal reflected by the plasmon-poiariton grating structure is defined 
as: 



r = - * 



(25) 



15 



where <p is obtained from the reflectance of the structure and c 0 and 1$ are the free space speed 
of light and wavelength, respectively. The group delay is obtained by applying finite- 
difference formulae of 0(h 4 ). This allows the numerical calculation of the derivative given 
discrete data points. 



DESIGN ARCHITECTURES 

It should be appreciated that the term ''cell" should not be interpreted to mean that the 
gratings disclosed herein necessarily are constructed from cells that are identical. In some 
gratings embodying the invention, the "cells" differ from one another. 

Uniform Periodic gratin g s 

The perturbation of the mode effective refractive index described above can be 
periodic and constant along the grating length. In this case, referring to Figure 3, the grating 
132 may be considered as a periodic perturbation comprising a series of cells 134, each 
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comprising two rectangular waveguides 100 and 100' having different lengths /, and l 2 and 
widths w, and w 2 , respectively. This cell 134 is then repeated N times to create the grating 
The dimensions of the waveguides in each cell 134, the spacing d, therebetween, the number 
of cells, and the spacing d 2 between cells are adjusted such that reflection occurs at a desired 
operating wavelength or over a desired operating bandwidth for an optical signal propagating 
along the grating axis, i.e. along the longitudinal axis of the array of cells 134. The period of 
the periodic grating is given by A„ = 1, + d t + 1 2 + d^. 

The center free space reflection wavelength of the grating is found using- 

A 0 = 2n ave \ 0 (26) 

where = w ' 0 ° /l0 ° * "' 0(1,/l00 ' i$ the average effective refractive index of the fundamental 

mode supported in each of the waveguide sections 100 and 100', n 100 and n 100 ' are the effective 
refractive indices of tbe.yj° mode supported by sections 100 and 100' respectively, 1 100 and 1 

are the lengths of the sections 100 and 100' respectively, and A<, is the fundamental period of 
the grating as defined above. For a fixed average effective refractive index, determined by the 
selection of the strip widths in the cell 134, the choice of the fundamental period will 
determine the reflection center wavelength. Design rules for the selection of appropriate 
lengths, // and l 2 , of the waveguides 100 and 100' respectively, include: 

. K 

(l.) Quarter-wave dimensions where the lengths, \ and 1 2 , are defined as l x = — and 



100' 



L = — - — where Xq is the center wavelength; 
4 "io<r 



(2.) Application of an algorithm, such as disclosed in reference [13], to find optimization 
parameters p opr and l opt for in the case of two lossy waveguide sections. This algorithm selects 
the lengths to maximize the intensity reflectance of the cell at the desired center wavelength. 
The lengths are then set to 1, = 0^1^ and J 2 = (1 - 3 opt )l opl . 

(3.) Application of a duty cycle. For a fixed fundamental period \, the cell subsection 
lengths, 1^ and 1 2 , are varied and design plots of Intensity reflectance with respect to grating 
length generated to provide design guidelines for dimension selection. The duty cycle is 
defined as the fraction of the metal strip length of waveguide 100 to the period length. The 
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duty cycle values of interest are limited by the minimum feature size of the targeted fabrication 
process. For illustration, consider a uniform periodic grating composed of a series of cells 
comprising a 8 ^m wide Au strip, 20 nm thick, and a gap filled by the background dielectric, 
Si0 2 . The fundamental period was set to 0.536 jim. Figures 6(a) and 6(b) show 2-D and 3-D 
5 plots respectively of the variation of the intensity reflectance with grating length and duty 

cycle. From these plots, it is clear that a tradeoff exists between grating length and duty cycle 
for maximum reflectance. The center wavelength can be found using equation (26). 

In the literature, the quarter-wave dimension and the 50% duty cycle cases are the 
most commonly used cell dimension definitions, 
lap The length L of the grating or N 9 the number of periods, is another design constraint, 

iph Figure 7 shows the variation of the intensity reflectance with N 9 the number of periods of the 
fl grating for a grating of the kind shown in Figure 4(b). The grating length and the number of 
f\ periods are related through the fundamental period as L = NA^. The performance shown in 
$ Figure 7 is for periodic uniform plasmon-polariton gratings with 50% duty cycle and period A 0 
l|1 = 0.536^im. The spacings d x and d^ in Figure 3 are set to zero. The width of waveguide 100 
111 W as varied from an 8 \xm wide Au strip to a 1 jim wide Au strip in increments of 1 ^m, 
O waveguide 100' is taken to be a gap filled by the background dielectric, SiO^ . The results 

show that an optimal length exists for each configuration, after which no additional reflectance 
is gained. From Figure 7, when the width of section 100 is set to an 8 \im wide strip, the 
20 optimal number of periods is N a 5000, which corresponds to a grating length of 2.68 mm. 
The existence of this optimal grating length is due to the losses that are present in the metal 
strip. 

The bandwidth of a grating AX is controlled by the width of the band-gap of the 
structure. In the case of a uniform periodic structure, the bandwidth is controlled by the 
25 effective refractive indices associated with the waveguide sections 100 and 100' of the cell 134 
and is defined as 

[f V 2 1 



JUL-26-2001 19=59 



ADAMS CASSAN MACLEAN 



1 613 828 0024 P. 28 



-19- 

where \ is the center wavelength, n avt is the average effective refractive index and An = n 100 - 
n m .\ is the modulation depth. Manipulation of the bandwidth thus requires a modification of 
the average effective refractive index of a cell. Techniques for the manipulation of the 
bandwidth will be described later in this disclosure. 

5 

Non-Periodic or chirped gratings 

The dimensions of the elements 100, 100' in each ceil 134 can also be made to change 
along the direction of the structure in order to implement a non-periodic perturbation of the 
mode effective refractive index over the grating length. The non-periodic perturbation can be 
a linear, quadratic or other function of position along the grating length. The grating is built 
ji from a cascade of N cells 134, were any of the parameters w h w 2> d h l„d 2 ,l 2 can vary 
£f between each adjacent ceil This creates a non-periodic or chirped grating design. The choice 
f\ of the widths in 134 will determine the depth of modulation. The lengths and spacing of the 
S cells will determine the center wavelength from equation (26). The lengths can be determined 
p using the procedure defined above, 

Pj Chirped grating designs are inherently non-reciprocal due to the presence of losses in 

rU the waveguide sections. Since the reflection of longer wavelengths requires longer cell 
n lengths, the long-wavelength reflecting sections are slightly lossier than the short-wavelength 
^ reflecting sections. Also, if the long- wavelength reflecting section is placed near the end of the 
20 grating, the propagation losses seen by the long-wavelength signals as they propagate through 
the entire structure result in a greater attenuation, compared to the short-wavelength signals, 
which are reflected near the input. The ordering of the cell cascade becomes important to 
equalize the attenuation over all wavelengths, 

25 Step-chirped gratings 

A step-chirped design is composed of the cascade connection of M different grating 
sections, each section being comprised of a uniform grating. The period of the uniform section 
varies from section to section to implement a desired chirping function. The utility of this 
design is to increase the reflectance created at each center wavelength defined by each uniform 

30 section. 



Mirror Interleave of non-periodic or chirped gratings 

The non-reciprocal nature of chirped designs can be mitigated via interleaving. A 
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mirror interleave is the process by which two gratings that are transposed longitudinally with 
respect to each other, Grating A and Grating B, each composed of cells of the form 1 34, 
referred hereafter as 134 A and 134B for Grating A and Grating B respectively, are combined 
in an alternating cell by cell approach. That is, a new grating is created by taking the first cell 
5 134A and cascading it with the first cell 134B, followed by the second cell 134A, and then the 
second cell 134B and so forth until all cells have been combined. The result is a quasi- 
reciprocal non-periodic grating design with improved spectral characteristics over the 
traditional chirped grating definitions in the presence of waveguide losses. See the preferred 
embodiments below for a detailed example. 

I 

£7 Interleaved gratings 

y*i The concept of cell by eel) interleaving can be applied to any of the design 

fp architectures disclosed herein. The interleaving is done as disclosed in the section captioned 
^ "Mirror interleave of non-periodic or chirp gratings", but no mirror or transposition is applied 
S to the gratings. This cell by cell interleave creates a new grating with a super-cell composed 
ry of the concatenation of one cell 134 of each original grating. The interleaving concept is not 
^ limited to the combination of two gratings; M gratings can be interleaved cell by cell in the 
M manner outlined. The resulting super-cell is composed of a minimum of 4 sections. The 
average effective refractive index of the super-cell is approximately equal to the weighted 
20 average of the average effective refractive indices of the individual cells or the weighted 
average of the effective refractive indices of the waveguide sections. 

Interleaving allows the positioning of the center reflection wavelength of the resulting 
grating in the spectrum, since a new average effective refractive index and fundamental period 
are obtained, and through equation (26), \ is subsequently modified. A shifting of the 
25 reflection spectrum is achieved by interleaving designs with different center wavelengths (that 
is, different dimensions /, and l 2 ). For example, using a 1 nm resolution on the cell dimensions, 
a grating was successfully centered at 1550 nm using super-cells each having 14-slices. 

The bandwidth of an interleave can be modified through the selection of the widths w 
composing the super-cell. Let us consider a super-cell of 4 slices (interleave of 2 gratings). 
30 For maximum reflectance, the modulation depth must be maximized between two slices. From 
the curves presented in Figure 7, maximum reflectance is obtained for the cell 134 
construction composed of a 8 \im wide Au strip 20 nm thick and a gap, where d t =d 2 = Q and 
50% duty cycle. The parameters, 4 / and w, of the remaining 2 slices are selected to fulfill a 
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center wavelength and bandwidth requirement, as the choice of the widths, and w 4) and 
lengths, /j and U referring to equations (26) and (27) affects the center wavelength and 
bandwidth via the average effective refractive index and grating period definitions. If the 
lengths are constant throughout the four slices, /, = l % = l 3 = /<> from equation (26), selecting 
5 (w 3i m>J<(w }7 will introduce a displacement in the center wavelength. 

Higher Order Architectures 

The device architectures described above were detailed considering only first order 
designs. The methodologies are not limited, however, to the design of first order 
fS embodiments. The center wavelength of higher order embodiments is found from: 

jS A„ = ^ (28) 

y § 

fj! where Mis an integer, greater than or equal to 1, that indicates the order of the grating period, 
T A 0 . For the same center wavelength and average effective refractive index, a higher order 
^ design is marked by a larger fundamental period since N> L Using equation (28) with any of 
PS the above disclosed design architectures, second, third or higher order gratings can easily be 
f4 designed into this plasmon-polariton grating technology. 

— The primary spectral characteristics of interest are the reflectance, the sidelobe 

suppression level, the bandwidth, the center wavelength and the group delay. The center 
wavelength and reflectance are manipulated according to the design schemes described above. 

20 Other design rules, detailed herein, permit further control of the bandwidth and sidelobe 
suppression. 



Suppression of sidelobe levels 

The traditional method for sidelobe suppression is apodization of the structure. The 

25 goal of apodization is to slowly vary the depth of modulation along the grating length, 

gradually increasing the coupling from the input and output of the grating, while keeping the 
average effective refractive index constant. This ensures that the center wavelength remains 
unchanged. Apodization functions found in the an can be applied to the invention to modify 
the bandwidth and sidelobe levels. 

30 An apodization is obtained, referring to Figure 3 grating 132, by varying the widths w 7 

and w 2 of the waveguides 100 and 100' to describe the slowly varying modulation depth. For 
example, w f would slowly increase while w 2 would slowly decrease from a nominal width 
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along the grating length. This creates a width contour or envelope along the grating length. 
In this scheme, the lengths of the waveguides are unchanged and hence the period is constant. 
This procedure requires that the cell 134 be constructed from two waveguide sections, that is, 
a gap cannot be used for either 100 or 100'. The apodized section created is then lossy. As 
5 these sections yield little reflection, they act as an absorber or attenuator at the input and 
output ends of the grating, decreasing the maximum reflectance possible. An alternative 
would be to use the background material for one partition of the cell 134. As the background 
is a lossless dielectric, this will reduce the absorption of the apodized sections and increase the 
maximum reflectance possible from the grating. This procedure requires fixing n m . to the 
ife background refractive index while n m is varied to implement the variation of the modulation 
% depth. The offset in the center wavelength perturbation is then compensated by slowly 
H* chirping the lengths /, and l 2 of the apodized section at a rate comparable to the modulation 
hj depth to ensure that equation (26) remains constant. This procedure can also be applied in the 
P ;i case of a two waveguide cell definition. In this case the waveguide with the smallest width 
1« would be held constant while the other waveguide width varies and a chirp is applied to both 
y dimensions i.e. the width and length. 

The apodization of gratings embodying this invention is dependent on possessing the 
O geometrical dispersion characteristics for a large number of waveguide widths with small 
incremental differences. It is impractical to find the modal solutions of a large number of 
20 waveguide widths so a cubic-spline interpolant of the geometrical dispersion characteristics of 
the symmetrical waveguide, as shown in Figure 2, is constructed. The interpolant of the 
normalized phase constant for a Au strip thickness of 20 nm is plotted in Figure 10(a), and the 
interpolant of the normalized attenuation constant is plotted in Figure 10(b). The diamond 
data points refer to the modal solution values; the curve is the result of the spline interpolation. 
25 This provides the required geometrical dispersion information for any desired width. 

Bandwidth 

The bandwidth of a grating is directly related to the band-gap of the periodic structure 
through the refractive indices of the sections, as shown in equation (27). In order to modify 
30 the bandwidth, a change in the average effective refractive index, n^, and/or the depth of 

modulation, An, must occur. Both parameters are dependent on the refractive indices of the 
cells, but the average effective refractive index is also dependent on the dimensions of the cell 
as it is a weighted average. The bandwidth can be reduced by increasing the modulation depth 
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and/or decreasing the average effective index, and vice-versa. This effect can be achieved by 
apodization if an appropriate trade-off between absorption and reduced reflectance is found. 
An easier scheme for bandwidth control is interleaving, This concept has already been 
disclosed in the previous subsections. It can be extended to any grating scheme to create a 
5 super-cell where the nominal bandwidth can increase or decrease according to the super-cell 
construction since n^ is perturbed, 

Center Wavelength 

The center wavelength of a grating design can be manipulated, according to equation 
M (26), by a variation in the average effective refractive index or the grating period. For a given 
p design, the choice of the widths and lengths of the waveguides 100 and 100' will determine the 
f|l center wavelength V To shift the location of \> using cells comprised of equal length 
S sections, an interleave with a second cell also comprised of equal length sections will cause a 
m perturbation in the average effective index while keeping the grating period constant, resulting 
lil in a shift in the location of V An interleave with a second cell of varying widths and lengths, 
f s [ will cause a perturbation in both the average effective index and the grating period, resulting in 
2? a more substantial tuning of V 

w In all of the above design methodologies, it should be noted that the waveguides 100, 

100* in each cell need not be rectangular, but a variety of other shapes could be used. For 

20 example, Figure 8 illustrates a portion, specifically two cells 138, of an alternative periodic 

structure 136 in which each cell 138 comprises two of the trapezoidal waveguide sections 106, 
106 ! with their wider edges opposed. It should be noted that these periodic structures, shown 
in Figures 3 and 8, are merely examples and not intended to provide an exhaustive detailing of 
all possibilities; various other periodic structures could be formed from cells comprised of all 

25 sorts of different shapes and sizes of elements. 

Substantially square waveguide strips 

The above device architectures can also be implemented using a plasmon-polariton 
waveguide structure having a substantially square cross-section strip that is finite in width, as 
30 is disclosed in co-pending United States Provisional patent application (Agent's docket No. 
AP877) filed June 22, 2001 . The substantially square waveguide structure, shown in Figure 
9(a), can be used to construct a new series of plasmon-polariton gratings composed of a 
metal film of thickness t, width where t is substantially equal to w } and equivalent 
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permittivity e 2 , surrounded by a cladding or background comprising an infinite homogeneous 
dielectric of permittivity t v Suitable materials for the waveguide strip include, but are not 
limited to, gold, stiver, copper, aluminum, and highly - or p-doped GaAs, InP or Si, while 
suitable materials for the surrounding material include, but are not limited to, glass, quartz 3 
polymer and undoped or very lightly doped GaAs, InP or Si. Particularly suitable 
combinations of materials include Au for the strip and Si0 2 for the surrounding material. The 
basic structure of the gratings based on this substantially square waveguide is shown in Figure 
9(b), which shows a grating 132C comprising a series of cells 134C, each comprising two 
substantially square waveguides 100C and 100C having different lengths /, and l 2 and 
substantially equal widths w ; and w 2 , respectively. This cell 134C is then repeated N times, as 
shown, to create the grating 132C. The dimensions of the waveguides in each cell 134C, the 
spacing d } therebetween, the number of cells, and the spacing d 2 between cells are adjusted 
such that reflection occurs at a desired operating wavelength or over a desired operating 
bandwidth for an optical signal propagating along the grating axis, i.e. along the longitudinal 
axis of the array of cells 134C. A perturbation of the mode effective refractive index can be 
established in at least two fashions in the grating 132C shown in Figure 9(b). For illustration, 
let us consider a uniform periodic grating, where the cell 134C is composed of two 
substantially square waveguides of differing metals, one of Au and one of Ag for example, in 
order to establish a perturbation of the mode effective refractive index. Alternatively, each cell 
134C can be composed of a substantially square waveguide 100C and a gap, being filled by the 
background dielectric or another dielectric, instead of waveguide 100C, where the length of 
the gap is adjusted for a desired operating wavelength. 

Tunabilitv 

It should also be noted that adjusting means may be provided for modifying an optical 
response of any of the gratings disclosed herein. For example, voltages can be applied to 
some or all of the strips in order to establish charges on the strips of the cells, which would 
change their permittivity and thus vary the optical transfer function of the structures of the 
device architectures disclosed above. Where the grating structure is as shown in Figure 4(a), 
for example, at least one electrode may be disposed alongside the grating and extending along 
at least a part of its length, as shown in broken lines in Figure 4(a). The electrode would be 
connected to one terminal of a voltage source and a second terminal of the voltage source 
would be connected to at least one of the strips. The distance between the electrode and the 
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grating would be great enough to ensure that optical coupling did not occur therebetween, 
while small enough to ensure sufficient field strength to effect the required modification of the 
response. If the dielectric material surrounding the strip were electro-optic, then the applied 
voltages would also change the permittivity of the dielectric, which also contributes to 
changing of the optical transfer function for a given design architecture. Where the grating 
has the structure shown in Figure 4(b), the adjusting means may comprise a voltage source V 
having one terminal connected to alternate ones of the waveguide strips and its other terminal 
connected to the intervening waveguide strips, as shown in broken lines in Figure 4(b). For 
further information about techniques for tuning the waveguide structures 100 used in 
embodiments of the present invention, the reader is directed to the afore-mentioned 
International patent application number PCT/CA00/0152S, 

Photonic bandgap structures can be created by placing two-dimensional arrays of cells 
(comprised of strips of various shapes and sizes) over numerous planes separated by dielectric 
material. The size and shape of the strips are determined such that stop bands in the optical 
spectrum appear at desired spectral locations, 

EXAMPLES OF PREFERRED EMBODIMENTS 

Unless otherwise stated, when structure dimensions are mentioned from this point 
onward, they refer to the Au/Si0 2 material combination, where the thickness of the Au strip is 
20 nm and the operating optical free-space wavelength is 1550 nm. The examples are of first 
order embodiments only. Dimensions and specification for higher order embodiments can be 
obtained by applying equation (28) to the same configurations disclosed. Similar dimensions 
are needed for most material combinations. 

The plasmon-polariton field may be excited by optical radiation coupled to the strip in 
an end-fire manner from a fiber butt-coupled to the input of the waveguide. The output of the 
waveguide can also be butt-coupled to a fiber, 

Example 1 : Uniform Grating 

Figure 1 1 illustrates an example of a uniform periodic grating embodying this 
invention. Grating 132D comprising a series of cells 134D composed of two waveguide 
sections, 100D and 100D', that take the form of an 8 \xm wide strip 100D and a gap 
lOOD'filled by the background dielectric, respectively, in this instance. The lengths of both 
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sections are selected for a 35% duty cycle and a 1 nm feature resolution where the period A of 
the grating is 0,536 \im. Figure 12(a), 12(b) and 12(c) illustrate the spectral performance of 
this embodiment for a device length of 4.288 mm. Figures 12(a) and 12(b) display the 
reflectance, linear and logarithmic, of the device, respectively. The uniform period 
embodiment yields a maximum intensity reflectance of 96.4%, with a 0.7 nm approximate 
bandwidth for a center wavelength of 1548,38 nm. Figure 12(c) displays the group delay in 
picoseconds of the uniform grating 132D over the useful bandwidth of the grating. As shown 
in Figure 12(c), the group delay is relatively flat over the useful reflection bandwidth of the 
grating 132D. 

Example 2 : Uniform Apodized Grating 

Figures 13(a) and 13(b) illustrate an example of an apodized uniformly periodic grating 
132E embodying this invention. It is composed of a series of ceils 134E each composed of 
two waveguide sections, 100E and 100£\ As shown in Figure 13(a), the contour of the 
© grating follows a predefined apodization function which varies the widths of the waveguide 
|ij sections 100E and 100E' as compared to a unapodized cell 134E shown in Figure 13(b), that 
= 5 takes the form of an 8 ^m wide strip and a 1 \xm wide strip, respectively. As an example, the 
M apodization function used is a sine profile given by 

/(*)=sm[yj (6) 

20 where L is the length of the grating and f[x) takes a value between 0 and 1 . The goal of 
apodization is to maintain the average effective refractive index constant even though the 
modulation depth becomes a slowly varying function of grating length. The sine profile 
perturbs the effective refractive index of the associated waveguides 100 and 100' according to 
equations (29) and (30); 

25 « 10 oW=^ + V/W (29) 

%W a " w -V/W ( 3 °) 

where h t = \n ave - n { | is the peak refractive index modulation, n flVe is the average effective 

refractive index and f(x) is the apodization profile. Once the required effective refractive index 
in each slice is known, the associated waveguide widths are determined to construct the 
30 grating. It should be appreciated that the cells 134E near the middle of the grating 132E will 
have substantially the same dimensions as the unapodized cell. of Figure 13(b). 
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Figure 14(a) and 14(b) show quantitatively the variation in the width and the effective 
refractive index, respectively, of the apodized-uniform periodic grating over the length of the 
grating 132E. The variation in the widths, w, and w 2 , of the waveguide sections 100E and 
100E\ imposed by the apodization profile is clearly shown in Figure 14(a) and the 

5 corresponding variation in effective refractive index in Figure 14(b). Figures 15(a), 15(b) and 
15 (c) illustrate the spectral performance of the grating of length 3 216 mm and period 0.536 
]Lim. Figures 15(a) and 15(b) display the reflectance, linear and logarithmic, of the grating, 
respectively. The present embodiment yields a maximum intensity reflectance of 80.5%, with 

Pi a 0,75 nm approximate bandwidth for a center wavelength of 1548.56 nm. A 22. 1 dB 
lfiM suppression in the sidelobes is achieved due to the sine-apodization profile. Figure 1 5(c) 

N= displays the group delay in picoseconds of the uniform apodized grating 132E over the useful 

y bandwidth of the grating. As shown in Figure 15(c), the group delay is almost flat over the 

J; : useful passband of the grating. 

* It should be noted that the invention is not limited to the use of sinusoidal apodization 

ill as in the above^described example but extends to various other forms of apodization profile, 

such as Gaussian, Blackman, hyperbolic tangent, positive hyperbolic tangent, and Raised- 
Cl Sinusoidal/Cosinusoidal, 

Example 3 ; Uniform period-wise interleaved grating 

20 Figure 16 illustrates an example of a uniformly periodic interleaved grating 132F 

embodying this invention. Grating 132F comprises a series of super cells 134F each 
composed of four waveguide sections, lOOF^ 100F 2 , 100F 3 and 100F 4 , that take the form of 
an 8 jim wide strip, a gap, a 1 \xm wide strip and a gap, respectively. The gaps are both filled 
by the background dielectric. The lengths L,, L 2 , L, and L 4 of the four waveguides sections 

25 100F b 100F a , 100F 3 and 100F^ , respectively, are adjusted such that L,~L 2 = 0.268 and 
L s = Lj = 0.269 ptm 9 for an overall period of 1 .074 |im. 

Figures 17(a), 17(b) and 17(c) illustrate the spectral performance of this embodiment 
for a grating length of 3.759 mm. Figures 17(a) and 17(b) display the reflectance, linear and 
logarithmic, respectively, of the grating. This uniform interleave embodiment yields a 

30 maximum intensity reflectance of 91 .8%, with a 0.4 nm approximate bandwidth for a center 
wavelength of 1551.15 nm, Figure 17(c) displays the group delay in picoseconds of the 
grating 132F over the useful bandwidth of the grating. As shown in Figure 17(c), the group 
delay is relatively flat over the useful reflection bandwidth of the grating 132F. 
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Exarnple 4 : Linearly chirped grating 

Figure 18 illustrates one example of a lineariy chirped grating 132G embodying this 
invention. Grating 132G comprises a series of cells 132G each composed of two waveguide 
sections, 100G and 100G', that take the form of an 8 \xm wide strip and a gap respectively, the 

5 gap being filled by the background dielectric. The lengths L x and L 2 of sections 100G and 
100G', respectively, decrease linearly between adjacent cells over the entire length of the 
grating 132G according to a chirp factor, C, The chirp factor C is equal to 0. Ixl0' 13 m . The 
lengths I, and L 2 of the first cell 132GJ are set to quarter-wave dimensions as detailed for a 
center wavelength of 1550 nm in the section captioned "Uniform Periodic gratings Figures 
l€f 19(a), 19(b) and 19(c) illustrate the spectral performance of this embodiment for a grating 

Jj length of 2.144 mm. Figures 19(a) and 19(b) display the reflectance, linear and logarithmic, of 

JsJ the grating, respectively. The embodiment yields a maximum intensity reflectance of 93%, 

III 

W with a 1 nm approximate bandwidth for a center wavelength of approximately 1 549.85 nm. 

y Figure 19(c) displays the group delay in picoseconds of the grating 132G over the useful 

iL bandwidth of the grating. As shown in Figure 1 9(c), the group delay is almost flat over the 

y useful reflection bandwidth of the grating. 



il Example 5 : Step Chirped grating 

Figure 20 illustrates one example of a step-chirp linear grating 132H embodying the 

20 present invention. Grating 132H comprises a series of cells 134H each composed of two 

waveguide sections, 100H and 1Q0H', taking the form of an 8 \im wide strip 100H and a gap 
100H', respectively. The gap being filled by the background dielectric. The lengths Li and L 2 
of sections 100H and 100H\ respectively, decrease linearly between adjacent periodic sections 
over the entire length of the grating 132H according to a chirp factor, C. The lengths L } and 

25 L 2 of the first cell 1 3 41^ are set to quarter- wave dimensions as detailed above, for a center 
wavelength of 1550 nm, in the section captioned "Uniform Periodic gratings". The grating 
132H is composed of 20 different uniform sections 135Hj, ... 135H 20 each comprising 200 
identical cells the lengths of which differ from those of cells in adjacent sections according to 
the chirp factor C =lxl0* 12 m. 

30 Figures 21(a), 21(b) and 21(c) illustrate the spectral performance of the grating 132H 

for a grating length of 2.1458 mm. Figures 21(a) and 21(b) display the reflectance, linear and 
logarithmic, respectively, of the grating 132H. The grating 132H yields a maximum intensity 
reflectance of 93%, with a 0,8 nm approximate bandwidth for a center wavelength of 1549.9 
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nm. Figure 21(c) displays the group delay in picoseconds of rhe grating 132H over the useful 
bandwidth of the grating. As shown in Figure 21(c), the group delay is relatively flat over the 
useful reflection bandwidth of the grating. 

5 Example 6 : Chirped period-wise interleaved gratings 

Figure 22 illustrates one example of a chirp-interleaved grating 1321 embodying the 
present invention. Grating 1321, shown in Figure 22, comprises a series of cells each 
composed of two waveguide sections, 1001 and 100I\ which take the form of an 8 wide 
strip and a gap, respectively , the gap being filled by the background dielectric. The series of 

10 cells is formed by a first set of cells A b A 3 , \ and a second set of cells A' s> A' 2 , A\. 

The two sets of cells are interleaved alternately, and are identical to each other, but transposed 
longitudinally one relative to the other. Thus, the grating 1321 shown in Figure 22 is the 
result of the period-by-period interleaving of the linearly-chirped gratings shown in Figure 
23(a) and 23(b). 

15 Figure 23(a) shows a linearly-chirped grating where the lengths, L f and L 2 , decrease 

linearly as the grating is traversed according to a chirp factor of C = 0. Ixl0' n m. The linearly- 
chirped grating of Figure 23(a) has 2000 periods. The lengths of its first cell are set to 
quarter-wave dimensions as detailed above for a center wavelength of 1550 ran. Figure 23(b) 
is the mirror image (flipping or transposing from left-to-right) of the grating of Figure 23(a). 

20 Figures 24(a), 24(b) and 24(c) illustrate the spectral performance of the grating 1321 for a 
grating length of 2.145 mm or 4000 periods, Figures 24(a) and 24(b) display the reflectance, 
linear and logarithmic, respectively, of the grating 1321 . The grating 1321 yields a maximum 
intensity reflectance of 92.9%, with an approximate bandwidth of 0.8 nm for a center 
wavelength of 1549.4 nm. Figure 24(c) displays the group delay in picoseconds of this grating 

25 1321 over the useful bandwidth of the grating, As shown in Figure 24(c), the group delay is 
almost flat over the useful reflection bandwidth of the grating. 



Summary of design methodology 

Figure 25 is a block diagram depicting the disclosed design sequence as applied to the 
design of all the preferred embodiments described hereinbefore. The design process begins 
with the definition of the waveguide structure to be analyzed as shown in the top left box of 
Figure 25. This information is used with our MoL formulation to solve for the dispersion 
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characteristics of the waveguides. As it is impractical to find the modal solutions of a large 
number of waveguide widths, a cubic-spline interpolant of the geometrical dispersion 
characteristics of the waveguide is constructed. These are the building blocks for the gratings. 
Starting in the top right box of Figure 25 the specifications for the grating and the design 

5 architecture of interest are selected. The 'cells' composing the grating are designed and the 
equivalent stack model of the grating is constructed using the results from the MoL 
formulation, A TMM approach is then used to estimate the spectral performance of the 
grating using the equivalent stack model. An iterative process can also be applied, as shown, 
to optimize characteristics. 
10*1 It is envisaged that the space between two strips could be either the same material as 

the background material, or a different material. For example, the space could be formed by 

H 1 etching a channel and filling it with the different material. 

f ; It should be appreciated that Figures 22, 23(a) and 23(b) show only five cells in each 

f ] set for simplicity; in practice there would be many more cells. 

II Although the specific embodiments shown in Figures 16 and 22 interleave the same 

%j kind of sets of cells to form the grating, it is envisaged that two different sets of cells, for 
l'i example a periodic set and a chirped set, to form the grating. 

O It should also be appreciated that the interleaved gratings could be fabricated using an 

alternative to the plasmon-polariton waveguide structures disclosed herein. For example, they 
20 could be fabricated using planar dielectric waveguides. 

While the invention has been described in detail in the foregoing specification, it will be 
understood by those skilled in the art that variations may be made without departing from the 
spirit and scope of the invention, being limited only by the appended claims. 
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